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DEC 77 L R SUTTON wee 
UNCLASSIFIED RASA/SRL-14-77-08 


ENR 
MED 


6-78 


B 


IFC(MFEA,EQ,0) Go TN 32 
TKNC(KT1)SSMLA(MS)®YKM(I,MS)#CTRE(L) 
KTLSKTLONRIFC 

32 IFC(MFLAP.EG,0) GN TO 34 
TKNCKT1)B3SMLA(MS)#CTRASCLIWYKMCI MS) 
KTISKTIeNRIFC 

34 IFC(MLEL,E0.0) GU TO 36 
TKNCKT1)SSMLA(MS) #@CTAG(L) a VKM( I ,MS) 

36 KTSKTo{ 

40 IF(MFEA,EQ,.0) Gu TO §0 
KT2SKT#NRIFC_- 
IFCHCT.E9,09) GN TO 4A 
TKNCKT) meF ABI (L) 


KT2mk T24¢NRIFC 

42 IFCMFLAP.EQ,0) GN TO 44 
TKN(KT2) BeFAB3Z(L) 
KT23K T24¢NRIFC 

44Y IFC(MLEL,EQ.0) GO TN 4b 
TKN(KT2)2eF ABY(L) 

46 KYBKT+} 

50 IFCMFLAP.E9,0) GO TO 56 
KT1akT 


IFCMCT,E9,.0) GO TU S2 
TKN(KT) @@FLBIL(L) 
KTISKTLI4NRIFC 

S2 IF(MFE4,EQ0.0) GO TO §4 
TKN(KT$) @eFLB2(L) 
KTISKTI4eNRIFC 

SG KTiIBKTI*¢NRIFC 
IFCMLEL.EQ.0) GO TU 55 
TKN(KT1) seFLBU(L) 

55 KTsxtel 

So IF(MLEL.E0.0) GO TN 60 
IF(MCT,EQ,0) GO TU 57 
TKNC(KT) meF SB1(L) 
KTSKT#NRIEC 

57 IF C(MFEA,EQ.0) GO TO $6 
TKN(KT) mer SB2(L) 
KTSKTONRIFC 

56 IF CMPLAP.EG,0) GO TO 60 
TKNCKT) geF SB3(L) 

( 60 CONTINUE 

900 CUNTINUE 
RETURN 
END 


547 
586 


SUBRUUTINE SURBMFC(1,J) 

INTEGER P,90 

COMPLEX EXPN 

CUMPLEX#16 TKN(441) 

CUMMON/INTER/NSY,NASEC »NFSEC,NByNBP,MFLAP,MFEA,MCT, 
IMFLEX»MCONsMAER,MFUS,NBC,»NFLAP,NFEA,NCT,NCON,NFFB, 
2NAS,/NHC,NVI,NSP,MAXN,NES, MSC »NEGN,IPCT»NITsMER,NORM, 
3IREM,NEX»NPSpNSCH,IG, IF ,NPRL»NPRS,NPD/NSK,NCOLS,/NCSH, 
4NFP1eMXSMIoMXT2P 1» MXKQeMXCPL pMXCSBeMXCPM,MXCPK »MXSMB, 
SNEBC,NESAC,MF ASH, MXFABsNFUS/NRAD,NRIFC pMXQ,NEIFCe 
ONEISC,NELTC,MY¥TKN,NFF,MINPN,MAXPN, IBF, MODE 

CUMMON/TKNI /TKN 

CUMMON/NLEAD/MLEL pNLEL »MCTY 

IF (NB,EQ.1)9G0 TN 586 

NMKaJel 

DO 547 MS22,NB 

DO 547 Qai,NRBO 

OU S47 Pa{,NRBD 

KTSNEIFC#NEISC#(Q01) eNRIFCONFUSSMXT2P1 oP 

KT2aKT¢(MSed )w(NEITCONRBD) 

TKNCKT2) 8 TKNC(KT) #EXPON(NMK yMS) 

RETURN 

ENO 


384 


11 


SUBROUTINE SOLVE 

REAL#8 DTLG10,DFAC,FAC 

COMPLEX#16 ZERO, TWO,S8UM,S8WAP,OTPHAS,DPIVOT,DETSV 
COMPLEX#16 EPS(63) FORCE (63) ,08(63,63) 
COMMON/EPSA/EPS 

COMMON /REMA/DETSV 
COMMON/CDOETRM/OPIVOT,DTPHAS,DILGIO,IOET 
COMMUN/DTERM/DFAC 

FAC#20,00 

TWORBDCMPLX(2.00,0,00) 
ZEROBOCMPLX(0,00,0,00) 

REWIND | 

READCL)MXSMISNRIFCeNHC »NORM,IREMY,NEX 
NORDERSMXSMI eNRIFC 

OO 11 Lut,MxSmI 

IKOWF aL aNRIFC 

TROwWSsTROWFeNRIFCH1 

OO 14 Kaei,MxSwHI 

ICOLF aK wNRIFC 

ICOLS@ICOLFeNRIFC#1 
READ(1)((COBCI,J),Ta#IRUWS,IJRUMF) »JBICOLS,ICOLF) 
REWIND 4 

OO 12 Isi,NORDER 

IKSNURDER*1@! 

FURCE (IK) SEPS(IK) 

DU 200 Jui,45 

WRITE (6,300) (DBC(I,J),121,45) 
FURMAT(5xX,10012,4) 

ILETsa4 

TSBLKSNHCwNRIFC 

ICOLS] SB8LK NORM 

TROWSISBLKSI REM] 

NEXCOL@I SBLKeNEX 

CALL SwWAPS(DB,NORVER,ICOL, TRUM) 
SwWAPBFORCECIROW) 

FORCE CIROW) SFORCE C(NORDER) 

FORCE CNURDER) sSwaP 

NSNORDER®@] 

CALL ERRSET(208,500) 

CALL OCmMAT(COB,N,FORCE) 

CALL ERKSET(208,10) 

WRITE (6,102) 0TPHAS,DTLGIO,0PIVOT 


102 FURMAT(//2X%,'DTPHAS @ ',2000.16e/2X,'DTLG10 @ ', 


1D40,16,/2x-.'DOPIVNT & ',2040,16) 
IFACSDOTLGI0O/F ac 

DFACSIFACwFAC 

OTLGLOSDTLG10eDF AC 
OTPHASSDTPHASa#DCMPLX(19, DO eel TLG10,0,00) 
DETSVSNTPHAS#®OPIVUT 

SUMBZERD 


DU 2 I81—9N 
IFCIZEQ,NEXCOL) GU TN 2 
SUMSSUM+DB(NORDER, I) eFORCE (I) 
2 CONTINUE 
SUMSFORCE (NORDER) =SuUM 
TF CCOARSCOBCNORDER,NEXCUL)) eNE.0eD9)GU TN § 
WRITE (6,6)NURDER,NEXCOL 
6& FURMAT('O DBC',IS," ,',1S,!' ) Ts Zero’) 
8TOP 
5 CONTINUE 
SUMSSUM/DRK(CNORDER,NEXCOL) 
FORCE (NE XCOL) S(SUMOFORCECNEXCUL) )/TWO 
FURCE (NORDER) @FORCECICOL) 
FURCECICOL)SZERU 
DO $ I8S1,NORDER 
3% EPS( I) SF ORCEC1) 
RETURN 
END 


SUBROUTINE SWAPS(DB,NORDER,ICUL + IROW) 
CUMPLEX#16 0B(63,603),SWAP 
DO 1 Ia1,NORDER 
SwAPaDB(T,ICOL) 
DBCI,ICOL) #0B(1I,NORDER) 

1 DBCI+NORDER)S3WAP 
DU 2 Il#i,NORDER 
SwAPSDB(TROW,I) 
DBC IROW,1)SDB(NORDER,1) 

2 DK(NORDER,1)8SwAP 
RETURN 
END 


651 
650 


100 
102 


104 


SUBROUTINE DCMATCALN,Y) 

REAL ®8 ADET,AMAG,USIGN 

COMPLEK#16 ApVY,0DET,NPIVOAT, TK yX 
COMPLEX#16 DATJ,ANX,DONE,DYI» TEMP, DAKJ »DYK pDAKK ,DAIK 
1,UNE,DPHAS 

DIMENSTIIN I6HG(63)2A063,63)0°¥663),X(65) 
COMMUN /CDETRM/ DPIVOTsDPHAS, ADET, IDET 
NDINS63 

ADE TsO ,00 

DSIGN21,0D0 

DPKASBICMPL X01,000,0,09) 

NP{sn 

IFCIDET,E.0) GO TO 650 

NPISN¢l 

DO 651 Isi,NP} 

X(T) SACNPY,I1) 

CONTINUE 

DU 118 Ke1,N 

AMX 3 A(K,K) 

IMXBK 

VU 100 TaK,N 

IFC CDABSCACI,K)) @LE, CDABS(AMX)) GO TO 100 
AMX # ACI,K) 

IMXs] 

CONTINUE 

IF CIMK,FQ,K) Gli TN 104 

DO 1904 Jai,,NP4 

TEMPSA(K,J) 

ACK, J)SACIMX,J) 

ACIMX, J) a TEMP 

ICHG(K)eIMx 

TEMPaSY(K) 

Y(K)s YCIMX) 

YCIMX)2 TEMP 

DPHASme[)PHAS 

GO frO 108 

ICHG(K) sk 

CONTINUE 

DAKKSBA(CK,K) 

FORMAT(1X,15¢29400,16/) 

WRITE (6,1000) DAKK 

FORMAT(SX, 'DAKK',5X%,2020,10) 
AMAGBCDABS(DAKK) 

IF (AMAG,NE,0.00) GO Th 6 

WRITE(6,7) 

FORMATC'O MATRIX IN DCMAT I8 SINGULAR') 
sToP 

CONTINUE 

ADE TBADE Te¢0L0G10(AMAG) 
DPHASSDPHASKDAKK /AMAG 
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an 


110 


653 
652 


112 


114 
116 
116 


120 
122 


124 


DONE SDCMPLX(1,00,0,00) 
DAKK@DONE /DAKK 

DO 110 Jat,NPy 

ACK, J) B8A(KoJ) eDAKK 
ACK, K) BDAKK 
IF(IDET,EQ,.0) GO TO 652 
TKaXx(K) 

DU 653 JmK,NPY 
XCJBXC Je TK ea(KeJ) 
CONTINUE 

DYKmY(K) 

Y(K) @DYK wDAKK 

DO 114 Jmi,Nn 

IF (J .EQ.K) GO) TO 114 
DAIKSACI,K) 


CALL ROWSUM(NPY,NDIM, ACI ,1),ACKo1) pDAIK) 


DO 112 Jai,NnP} 

ACI, JIBACI oJ) mDAIKRA(KeJ) 
ACI,K)SDAIK 

OYIBY(T) 

OVKaY(K) 

Y¥CLT BDV LeDALK eDYK 
CONTINUE 

ACI ,K) meA(I,K) RDAKK 
ACK, K) SD AKK 

CONTINUE 

OO 122 Ke{,N 

LaNelek 

KI@ICHGE(L) 

IF (L,EQ.KI) GO TU 122 
OU 120 Imj,N 

TEMP @ ACT>L) 

ACI,L) w ACIekT) 
ACI,KI) @ TEMP 
CONTINUE 

IFCIDET.NE.0) DPIVOT#Xx(NPt) 
RETURN 

ENO 
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COMPLEX FUNCTION EXPON(L,MS) 
C CREATE EXP(I wl wPHIM) 

DIMENSION CS(U,6)-SN(4,6) 
COMMON/RNAM/CS,SN 
ILSIABS(L) 
IF(L) 160,15,17 

15 EXPONSCMPLX(1,.0,) 
GO TO 18 

16 AaCS(MS,IL) 
BRSNCMS,IL) 
EXPONECMPLX(A,@6) 
GO TO 18 

17 AgsCS(MS,L) 
BESN(MS,L) 
EXPONSCMPLX(A,B) 

18 CONTINUE 
RETURN 
END 


17 


21 


20 
30 


SUBROUTINE SwAa(I,J) 

INTEGER P,Q,98 

COMPLEX EXPON 

COMPLEX#16 CMSeCM1,CS810CS82 

COMPLEX#16 TKN(441) 

CUMPLEX#16 ZLN 

CUMPLEX#16 XNLQ 

COMPLEX#16 YKM(3,4) 
CUMMON/SPAR/AKCI (6), TAUC6),SMLA(6) eDMS(6),AK(U9,AC(CU), BUC), 
1CAPK,CAPC , 
COMMON/INTER/NSY,NBSEC,NFSEC,NB,NHP,MFLAP,MPEA,MCT, 
IMFLEX»MCON,MAER,MPUS,NBC,NFLAP,NFEA,NCT,NCON,NFFS, 
2NAS,NHC,NVI,NSP,MAXN,NES,MSC eNEGN,IPCTe¢NITSMER,NURM, 
SIREM,NEX»NPSsNSCHS IG, IF »NPRL »NPRS,NPD,NSK,NCOLS,NCSB, 
GNFP 1g MXSMI gMXT2P 4 oMXKOOMXCPL »MKCSHMXCPM,MXCPK »MXSMB, 
SNEBC,NESBC,MFASHsMXFAByNFUS/NRBD,NRIFC pMXQ,NEIFC, 
ONEISC,NEITC,MXTKN»NFF,MINPN,MAXPN, IBF pMONDE 
CUMMON/NLEAD/MLEL eNLEL oMCTY 

COMMON/ROTF /0M1,0M2,0MT 

COMMUN/FREF/CMS,CM1,CS1,CS82 
COMMON/SWASH/SWGJ,SWEI, SWM,SWR 

COMMUN/STKNI/STKN 

KSNLaSIeNFPI 

IFCI,NE.J) GU TO 30 

DO 17 Lai,MxTePi 

LSZLeMAXNel 

LLe(Lel) eNRIFCOoL 

TKN(LLIMBZLN(LS,T) 

DO 20 4S$a1,NB8 

DU 20 Lal,MxTaPy 

LSSLeMAXKNeal 
CFOL SL, OFOMS(MS) w(1L o(L Sel Sel) /(1, oL Sal SaSwGJ/SwE?))/SWR 
LLe( Le] ) eNRIFCOMXT2P14¢(MSe{ ) *NRBDONCOLSENFUS 

IF CMPLEX.EQ,0) GN TO ai 

IF (MCON,EQ.0) GO TO 20 

LLSLL¢3 

IFCMCTY,GT.0) LLSLL?3 

TKN(LL) SPE XPON(L S98) #CFOL 

GO TO 20 

LLSLLeMCT 

IF (MCT,.EQ,0) GO TO 20 

CSISCMSeCM1i ak SMLeO0M} 
YKMCI,MS)@SMLA(MS)&(1¢CS{ x TAU(MS§)) 

TKNCLL I @YKM(I,MS)XEXPON(LS,MS) @CFOL 

CONTINUE 

GO TO §0 

IMJaleJ 

DO 18 Lwi,MxTePi 

DO 18 Qm1,MxT2P4 

IF(L.E@,9) GO TO 18 


LMOsLeQ) 
IFCIMJ,NE.LMQ) GN TO 18 
LSZLeMAXNel 
ASZQeMAXNel 
LLS(Le1)eNRIFCSQ 
TKN(LLIAXNLOCI>LS,Q8) 
156 CONTINUE 
50 RETURN 
END 


SUBROUTINE SWB(I,J) 

INTEGER P,Q,08 

REAL®6 CX(75) 

COMPLEX EXPON . 

COMPLEXw16 ULN,S(216) 

COMPLEX#16 TKN(441) 

COMPLEX#16 B(648),SM_L8(106),8MLC(6108),SMLD0(108) 

COMPLEXwi6 CTA(S4),FAB(SU),FLB(54) 

COMPLEX#16 CTB1(9) -CTB2(9),CTB3(9),FABL (9), FAB3(9) pFLB1(9) FLB2(9) 

COMPLEXwi6 SMLE(108),FS8(54) 

COMPLEXw16 SMLF(108),SMLG(108) 

CUMPLEX#16 FSB1(9),F882(9) ,F$83(09),FAB4(9),FLB4(9) CTB4(9) 

CUMMON/BTS/By¢ SMLB, SMLC + SMLO,CTB,FAB,FLB,CTB1 +CTR2,CTB3eFABI, 
1FA83,FL61,FLB2,SMLE,FS8,FSB81,FS$82,F883,FAB4,FLBG,CT84,SMLF,SMLG 

CUMMON/SWASH/SWGJ» SWEI » SWM,SWR 

COMMON/INTER/NSY,NBSEC»NFSEC »NByNBP,MPLAP,MFEA,MCT, 
IMPLEX,MCUN,MAER,MFUS,NBC,NFLAP,NFEA,NCT,NCON,NFFO, 
2NAS,NHC,NVI,NSP,MAXN,NES,MSC pNEGN,IPCT»NITsMER,NORM, 
ZIREM,/NEX»NPS,NSCHe IG, IF,NPRL»NPRS,NPD,NSK,NCOLS,NCS88, 

UNFP 1» MXSMI pMXT2P1 eo MXKQoMXCPL eMXC8BYMXCPM,MXCPK MX SMB, 
SNEBC »NESBC»MFASB,MXF AB, NFUS,NRBD oNRIFC oMXQeNEIFCe 
ONEISC,NEITTC,MX TKN »NFF,MINPN,MAXPN, IBF, MODE , 

COMMON/SPAR/AKCI (6), TAUCO) » SMLA(6) ,DMS(6),AK (49, ACCU) BUCK), 
1CAPK ,CAPC 

CUMMON/CFLEX/CX 

CUMMOUN/TKNI/TKN 

COMMUN/NLEAD/MLEL oNLELoMCTY 

CUMMUN/881/8 

NONUS] 

IMJaleJ 

JMNCB( Jeol) ei2eNCOLS 

IF CMFLEX.EQ,0,AND,IMJ,NE.0) GO TO 23 

IF (MFLEX.EQ,0) GO TO 3 

IFCMCTY,€0.0) GO TO 3 

IFCIMJ,NE.0) GO TO 23 
3 IFCIMJ,NE,0) NONOSO 

IF (NBP,E9,0) GO TU 13 

KSMLOIeNFP] 

NPMKB@NP Sek SML 

DO 14 Qei,MxT2Py 

QSBQ@MAXNe} 

NMKGSNPMK @QS§ 

IFC(NMKO,EQ.0) GO TO § 

RNMKQ@1, OeNMKG 

RF ABRNMKQ/NBP 

NF ABSNMKQ/NBP 

DIF MRF Aw], OwNFA 

IFCUIF,GE.0,0) GO TO 2 

DIFseDIF 
2 IFCOIF,GT,,.05) GO TO 14 
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14 


13 


24 


25 


22 


21 
23 


DFGR#1L,OFDMS( 1) 01,4 CUSHOS—1,)/01, +08xOSeSWGJ/SWET))/SWR 
IF CMPLEX.E9,.9) GO TO 7 

IF (CMCON,EQ,.0) GO TO 14 

IFCMCTY,GT.0) GO TO 6 

OU & I[O81,NCOLS 
LLENEIFC#NESCe( LOL) eNRIF CHO 

LeC( Jel) #NCOLSaMXSMIe(ITeL I eNCOLS+IQ 
TKN(LL) SeNBPROFQR*FLACL) 
LLBNEITFCe#NELSC#NCOL SeNRIFCHO 

TKNOLL) SB ¢NBPeDFOR#C X03) wNOND 
TKN(LL#NRIFC) s#NBPwOFQR&CX (8) *NOND 
TKN(LLEONRIFC#2) seNBPwDFQR*CX (12) *NONOD 

GU TO {4 
LLENEITFCe#NEISCHeNRIFCaeNCOL S4+S#NRIFC HQ 
TKNC(LL) seNBPeDFOR 

GO YO 14 

LLENETFC#NEISCeNRIFCHNCOLSOQ 

IF CMCT.EQ,0) GO TO 14 

TKN(LL) aeNBPRDFOR/SMLAC1) 

CUNTINUE 

GO TO 23 

DU 21 MSe1,NB 

DU 21 Pal,MxT2P4y 

QSSPeMAXNel 

DFQARELS+OMS(MS) e(14(O8Se08—1)/(1,+0SnOS*SWNGI/SWEI))/SWR 
IF CMFLEX,E9,0) GO TO 22 

IF(MCON,EQ,.0) GO TO ai 

IPCMCTY,G6T.0) GO TO a5 

DU 24 ITQay,NCOLS 
LLENELFCeNETSCH(MSe1 RNETTCo (CI Ge] )enNRIFCoP 
LeCJel) eNCOLSeMXSMI4( Jw} )eNCULS+IQ 
TKN(LL) SeE XPON (@QS,MS8) *OFQReFLBE(L) 
LLENETFC#NEISC#(MSa1) eNETTC+NCOLS#NRIFC OP 
TKN(LL) @¢EXPON(@QS9MS8) wDFOR*CX(3) eNONO 
LLSLL#NRIFC 

TKN(LL) RFE XPON(@9S,M8) DF URC X (6) eNONO 
LLELLENRIFC 

TKN(LL) @wEXPON (09S 5M8) *DFQRRCX (12) *NOND 
GU TO 21 

LLENEIFC#NEISCH#(MSel )eNETTCONRIFCHNCOLS +S aNRIFC OP 
TKN(LL) mw E XPON(@Q3.M8) ®UFOR 

GU TO 21 

LLSNEIFC#NEISCO(MSa] DeNETTCeNRIFCHNCOL SeP 
IFCMCT.EQ.9) GO TO 24 

TKN(LL) BOE XPON( @NS,M9) RDFOR/SMLA(MS) 
CONTINUE 

CONTINUE 

IF (MFUS,EU,0) GO TU S52 

JMI sel 

DU 50 IPPai,MxT2Pi1 
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LSS] PPeMAXNeil 

IFCUMI.NE.LS) GO TO §0 

DU S1 IuQei1,NCOLS 

LLPSIMNCO(1QQe) Jaleo} 

LLENEIFCo(TQ0e1) aNRIFCOI PP 
31 TKN(LLISULN(LSsKSML) @S(LLP) 
50 CUNTINUE 
Se CUNTINUE 

RETURN 

END 
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10 


i} 


COMPLEX FUNCTION ZLNw16(LS8,1) 

REAL*8 C3,Cu 

COMPLEX EXxCHI 

CUMPLEX«16 CLNJ,C1.C2,C5,C6,C7,C8,C9,C10 

COMPLEX#16 CMSeCM1,CS812CS2,VNe VNI,VN2 
CUMMON/SWASH/SWG6J,SWEI,SWM,SHR 
COMMON/INTER/NSY,NBSEC/NF SEC eNB/NBP,MFLAP,MFEASMCT, 
IMFLEX»MCON,MAER,MFUS,NBC,NFLAPsNFEA,NCT»/NCON,NFFB, 
2NAS,NHC,NVI,NSP,MAXN,NES,MSC,NEGN,IPCTONIT»MER,NORM, 
ZIREM,NEXeNPSeNSCHP IG, IF ,NPRL e¢NPRS,NPD,NSK,NCOLS,NCSB, 
UNFPLgMXSM] -MXT2P1 yMXKG,MXCPL »MXCSB»,MXCPM,MXCPK MXSMG, 
SNEBC »NESAC »MFASB,MXFABe,NFUS,NRBD ,NRIFC pMXQONETFCe 
ONEISC »NELTC,MXTKNeNFF,MINPN,MAXPN, IBF p MONE 
CUMMON/SPAR/AKCI (6), TAU(H) » SMLAC HO) »DMS(6),AK(U),A4CCU),BU(4), 
1CAPK,CAPC 

COMMON/NLEAD/MLEL op NLEL oMCTY 

COMMON /FREF /€MS,CM1,C81,CS2 

COMMON/ROTF/0OM1,0M2,0MT 
COMMON/SPARL/AKT(4),ACT(U) pAKPC4), ACP C4) 

RsSwWR 

KSMLBIeNFP{ 

CS1aCMSaCM1 eK SML eM) 

CS2aCSixCsi 

CISSWMe(CS2eOMTaLSeCMieCStleL Set Same) 
C2eDCmMPLx(0,00,0,00) 

C9BDCMPLX(0,.00,0,00) 
CFLB1.+(LSeL Sel) /C14L SeLSuSwGJ/SWET) 

CPLRSCFL/SWR 

DO 10 JJsi,NES 

CLNJBAK (JJ) CSI WAC (CIS OCMILMLS#OM1nAC (JJ) 

Ceace+C ing 

COBCIOSCLNJ a (1 eB (JJ) eCFLR) 

C10aCe 

Cesc9 

IF(MSC,EQ,0) GO TO 11 

Cescio 

CSBCAPKSCAPCa(CS16CM1 wl Se0M1) 

Cos0CMPLx(0,00,0,09) 

C7S0CHPLX(0,0050,00) 

DO 6 JJa1,NES k 

COMAK(C JJ 4CSin ACCIJ MCh ML SwOMIWAC (JU) 
COBCOSCBWEXCHI(LSO,JJIaC1 eB CUI) WCFLR) 
CTSCTH+CBREXCHICOsLS, JJ aC leBUC JJ eCFLR) 
C28C9el74C6/(CS5+C2) 

TE CMAXN,EQ.3) GO TO 12 

KXS1@LSelS 

C382 ,00835, 141592654 ul Sa Sak XaKX 
CGBReRaRa (1, /SWCJPL Sel S/SwEI) 

C38C3/C4 

ZLNSCLOC2e¢DCMPLX(C3,0,00) 
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le 
13 


GO TU 13 

ZLN32C14C2 

CUNTINUE 

VNEDCMPLX(0,.00,0,00) 

VNISDCMPLX(0200,0,.00) 

VN@SUCMPLX(0.00,0.90) 

DO 5 JJsi,NES 
VNSVNGAKT(JJI4+(CSLECML aL SaNMt )wACTC SJ) 
VNLBVNICAK (SII + (CSL eCM1 wLSwOMl )wAC (JJ) eB CJ eC ReBI (JJ) CFL) 
VN2BVNAGAKP (JJ) 4 (CS10CM{ wLSOM1) ACP (UJ) 
ZLN@ZLNG(CLSaLSeVNOVN{ ACFLE¢CFLACPL aVN2) J(ReR) 
RETURN 

END 
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COMPLEX FUNCTION XNLQwi6(1,L39e9S) 
INTEGER QS 
COMPLEX EXCHI 
COMPLEX#16 C81,CS2eCM1 ep CMS,XNoXN1 yp XNSAXNS 
COMPLEX#16 CS,WN,WN1,WN2 
COMMON/SWASH/SWGJ» SWEI» SWM, SWR 
COMMON/INTER/NSY,NBSECsNFSEC pNBpNBP,MFLAP,MFEA,MCT, 
IMFLEX,MCON,MAER,MFUS,NBC,NFLAP,NFEA,NCT,NCON,NFFE, 
2NAS,NHC,NVI,NSP,MAXN,NES,MSC,NEGN,IPCT»NIT»MER,NORM, 
ZIREM,NEX,NPS,NSCH, IG, IF ,-NPRL»NPRS,NPD,NSK,NCOLS,NCS8B, 
UNFP 1» MXSMI pMXT2OP1 > MXKO,MXCPL »pMXCSB,MXCPM,MXCPK ,MXSMB, 
SNEBC,NESBC,MFASK,MXFAB, NF US», NRBDNRIFC pMXQ,NEIFC? 
ONEISC,NEITC MX TKN, NFF,MINPN,MAXPN, IBF,» MODE 
CUMMON/SPAR/AKCI (6), TAU(6) »SMLA(6) »DMS(6),AK(4U),4C€4),B8304), 
1CAPK,CAPC 
CUMMON/NLEAD/MLEL oe NLEL po MCTY 
CUMMON/FREF/CMS,CM1,C81,CS2 
COMMUN/ROTE/OM1,0"2,0MT 
COMMON/SPARY/AKT(4),ACT(U),AKP(4),ACP(4) 
KSML@IeNFPIL 
CSISCMSeCM {eK SML eUM1 
IFC QS,EQ.LS) GO TU 15 
XNLOBSDCMPLX(0,00,0.00) 
CFQ81, (9980801) /(1, 40S eOSeSWGJ/SWET) 
CFL EL +(LSwLSal)/(1FLS*LSaSwGJ/SWET) 
CFQRECFQ/SWR 
CFLRSCFL/SWR 
DO 10 JJa1,NES 
10 XNLUSXNLO+CAK( IJ) *(CS1°CM1eOSROM1) RAC(C US) WEXCHICLS,OSeJJ) & 
1¢1,e83 (JJ) aCF QR) 
IF (MSC.E%9.0) GO TO 16 
CSBCAPK+CAPCH(CS10CM1 eOSe0M1 ) 
XN1BDCMPLKX(0.00,0,00) 
XN2SDCMPLXK(0,.00,0,00) 
XN3SRBDCMPLX(0,00,0,00) 
DO 12 JJ#1,NES 
XNBSAK(JJI¢(CCSLSCM1eOSHmM1) wAC (IJ) 
XNIBXNLOXN 
KXN2RXN2¢XNWEXCHI( LS» Op JJ aC LeeBIC JJ) MCFLR) 
12 XNSSXNSOXNMEXCHICO,OS,JJI#(1.°80 (JJ) 2CFOR) 
XNLOSXNLQeXNSRKXNA/S(CS4XN1) 
GO TU 16 
15 XNLOSDCMPLX(0,00,0.00) 
16 CONTINUE 
RweSwR 
WN BDCMPLX(0,D0,0,00) 
WNISDCMPLX(0,00,0,00) 
AN2BD0CMPLX(0,09,0,00) 
00 S JJsa1,NES 
WNBWNOC(CAKT( JJ) *(CSLOCMI ROSHUM1)MACT( JJ) ) WEXCHI(LS,Q8,JJ) 
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Re er mee 


WNISMNLOCAK (JID (CCSLeCMLeQSHOMd WAC (JJ) BIC JI eC ReBICIIIRCRA)® 
1 EXCHI(LS,QS-JJ) 

5 WN2MWN24¢(aKP( JJ) +(CS10CM1 eQSHUMS ) ACP (JJ) ) REXCHICLS,OSeJJ) 
XNLUBSXNLUS COSRLSHWNeWNI RCRLECKQRCEL eWN2) /(RAR) 
RETURN 
END 


COMPLEX FUNCTION ULNe16(LS,KSML) 

CUMPLEX«16 UN, UNNZUNC,CS1,CS2,CMS,Cul1 

COMPLEX EXCHI 

CUMHON/SPAR/AKCI (6), TAU(6) » SMLA(6) »OMS(6),AK(U), ACCU), BUCY), 
{CAPK,/CAPC 

COMMONSROTFE/041,90M2,9MT 

CUMMUN/FREF /CMS,C41,081,CS2 

CUMMON/SHWASH/SWE J, SWEILZSHM,SWR 

CUMMON/INTER/NSY,NBSEC pNFSEC,NAeNBP,MPLAP,MFEA,MCT, 
IMFLEX,MCUN,MAER,MFUS,NBC »NFLAP,NFEA,NCT»NCON,NFFB, 
QNAS*NHE pNVIZNSP,NAXN,NES,MSC pNEGN,IPCToNTT»MER,NORM, 
SIREMPNEX»NPS*NSCHS IGS IF eNPRL»oNPRS,NPD,NSK,NCULS,NCSB, 
UNFPL gp MXSMTpMXT2P Lp MXKOMXCPL eMXCSEPMXCPM,MKCPK MXSMB, 
SNEBC,NESBC, FASB, MXF AB, NFUS,NRBL,NRIFC pMXQ,NEIFCS 
ONEISC,NETTC,MXTKN eg NEFF »MINPN,MAXPN, TBF ,MODE 

CUMMUN/NLEAD/HLEL @NLEL MC TY 

UN SDCMPLxX(0,00,0.00) 

UNNSDCMPLX(0,00,0,00) 

CFLSL,+(LSaL Sal) SCL eLSeLSaSnGJ/SWET) 

CFEFLRSCFL/SWR 

CSISCHSeCM1 eK SML *UM] 

Du S JJz1,NES 

UNSUN+ (AK (JJ 4¢(CSILECMIALSHOM1) RAC( JJ) MEXCHICOOLS,Jud® 
1C1,eHU (JJ) eCFLR) 
S UNNSUNN@(AK (JS) + (CSL eCMLeL Salmi wal (JJ)) 

IF(MSC,EQ,0) GO TU 6 

UNCSCAPKS(CS1eCMi AL Sah) RCAPC 

ULNSUN®UNC/(LINCeLINN) 
6 TFCHSC,E9, 9) ULNSUN 

RE TURN 

END 
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15 


16 


17 


18 


COMPLEX FUNCTION EXCHI (L,9.J) 
INTEGER Q 

CUMMON/RNAM] /CS8A(4,20),8NA(4,24) 
LOsLe@ 

IL@aIABS(LQ) 

IF (L0)16015,17 
EXCHISCMPLX(1,.0,0,0) 

GO TO 18 

AeCSA(J,ILQ) 

BsSNA(J,ILQ) 

EXCHISCMPLX(4A,eB) 

GO TO 18 

AsCSA(J,ILQ) 

BsSNA(J,ILQ) 

EXxCHIasCmMPLx( 4,8) 

CONTINUE 

RETURN 

END 
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SUBROUTINE ZTEGICI,J) 

INTEGER P,Q 

COMPLEX EXPON,EXPM1,@XPP1 

COMPLEXe16 Cm) 

COMPLEX#16 B8(648),8ML8(108),8MLC(108),8MLD(108) 

CUMPLEX#16 CTB(54),Fa8(S4),FLE(S4) 

COMPLEXwi16 CTB1(9) ,CTB2(9) CTE3(%) ,FAB1(9) ,FAB3(9) , FLO1(9) .FLB2(9) 
COMPLEX#16 $MLE(108),F88(54) 

COMPLEX#16 SMLF(108),SMLG(108) { 
COMPLEXa#16 FSB81(9).F$82(9) .F883(9),FAB4(9),PLO4(9),CTB4(9) 
COMPLEX#16 TKN(441) 

COMMON/BTS/8,SMLB,y SMLC,SMLO,CTS,FAB,FLE,CTBI »CTB2,CTBS,FABI ,FABS, 
1FL61,FLB82,SMLE,FS8,F881,FS82,FSB3,FAB4,FLB4,CTB4,SMLF,SMLG 
CUMMON/TKNI/STKN 
COMMON/INTER/NSY »NBSEC ,NFSEC,NB,NBP,MFLAP,MPEA,MCT, 
IMFLEX,MCON,MAER, MFUS,NBC,NFLAP,NFEA,NCT,NCON,NFFB, 
2NAS,NHC NVI ,NSP,MAXN,NES.MSC oNEGN, IPCTeNITSMER,NORM, 
SIREM,NEXsNPS,NSCH/IG,IF,NPRL »NPRS,NPD,NSK,NCOLS,NCSB, 
GNEP 1» MXSMI »MXT2P1 oMXKQ>MXCPL »MXCSB,MXCPM,MXCPK,MXSMB, 
SNEBC ,NESBC MF ASB MXFAB,NFUS,NRBO NRIFC pMXQANEIFC, 
ONEISC pNEITCeMXTKNGNFF,MINPN,MAXPN, IBF, MODE 
CUMMON/NLEAD/MLEL oNLEL pMCTY 

Cmis DCMPLX(0,00,3,00) 

NMKse Je! 

KBS NEITFCONEISCeMXT2P14¢NFUS 

KBBS KBeNCOL SaNRIFC 

NMKPISNMK eo} 

NMKMLaNMKel 

NPKSZeNPSel] oNFP{ 

NBM18(NBe] ) #NRBD 

LSMAB( Jol) mNESBCe(10}) eMxC88 

LLARS( Jol) wNEBCe(Iel])aMXCPM 

IF (MFLEX.EQ@.0) GN TO 10 

WRITE (6,900) 

FORMAT(/,9X,'GIMBALLED OR TEETERING ROTOR, MFLEX MUST EQUAL ZERO') 
GO TO 90 

IF (NBC,EQ,2)GO TN 13 

IF (NBP.EQ,0)GU TO 11 

IF(NBP,LE.2)G60 70 i2 

GO TO 16 

IF(NB,GT.2)G0 TO 16 

WRITE (6,901) 

FORMAT(/,9X,'GIMBALLED ROTOR MUST HAVE MORE THAN TWO BLADES!) 
GO TO 90 

IF (NBP.EQ,0)G0 TO 14 

IF (NBP,NE.2)GO T0 15 

GO TO 16 

IF(NB,EQ.2)60 T0 16 

WRITE (6,902) 

FORMAT(/,9X,"TEETERING ROTOR MUST HAVE TWO BLADES!) 


GO TU 90 
fo NES 
NAS] 
IF CNHC,EQ,1960 TN 17 
NE St 
NAS) 
17 Nise! 
N2Z6 
DO @2 NNw{,NE 
IF(NN,FQ.1)9G0O TO 18 
N131 
Nese 
18 DO 19 Qst,NCOLS 
LSSLLAR4+(Qel)el2e3 
L1i0=sL 347 
KKSKA4¢(Qel )eNRIFCON2 
19 TKN(KK) SNARBA(L3)ONIRCMIRB(L10) 
L3SLSMAe3 
L1i0SL 347 
KK SKAB4¢NQ2 
IF (MCT,E9,0)60 TO 20 
TKN( KK) SeNARSMLB(L3) eNI CMI SMLB(L10) 
KKSKK#NRIFC 
20 IFC(MFEA,EQ,0IG0 TU 21 
TKN(KK) seNARSMLCO(L3)@NIRCMI®SMLCE(L10) 
KKSKKONRIFC 
ei IF CMFLAP,EQ,0) GO TU 23 
TKNCKK) BeNARSMLD(L3) eNIRCMIRXSMLDO(L10) 
KKSAKONRIFEC 
23 IFCHLEL,EQ.0) GN TO 22 
TKN(KK) SeNARSMLE(L3)@NIWCMIRSMLE(LI0) 
22 CONTINUE 
IF (CNBP.NE,0)GN 10 40 
TF CNBC,£9,2)G60 TO 28 
OW 24 Qai,NRBN 
KTS3KB84(Qel)wNRIFCH2 
KKBKT44 
K2SK TONAM 
KOSKKeNBMY 
TKN(K2) B@TKNCKT) 
24 TKN(K6) Be TKN(KK) 
DO 25 MS$a2,NB 
EXPM1SE XPON(NMKM1,MS) 
EXPP1SEXPON(NMKP1,MS) 
MSHIF T8(MSe1) wNEITC4(MSe2) @NRBD 
DO 25 Qsi,NRBD 
KTSKB4(Qej)aNRIFC+?2 
KKBK TG 
K28KTOMSHIFT 
KOBKKOMSHIFT 
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es 
2a 


30 


31 


$2 


34 
33 


35 


38 


40 


42 


K22aK 2¢NRB0 

KO6SKO4eNRAED 

TKN(K2) weTKN(KT) wEXPR) 
TKN(K6) BeTKN(KK) wEXPM 
TKN(K22) Me TKN(K2) 
TKN(K66) Ze TKN(K6) 

GO TO 100 

00 30 Qst,NCOLS 
LIISLLAR+(Qel)wlarti 
KKSKB¢(Qe1 )eNRIFC#NRADS6 


TKN(KK) SB(L41) 
LilsLSmMaei} 
KKSKBBeNRAD 6 
IFCMCT.E9,0)G0 TO 31 
TKN( KK) @@SMLB(L11) 

KK SKKENRTFC 

IF (MFEA,EQ.0)G0 TO 32 
TKN( KK) me SMLO(L11) 
KKSKKE#NRIFC 

IF (CMFLAP,EQ,.0) GO TO 34 
TKNC KK) @@SMLO(L11) 

KK SKK ONRIFC 


IF CMLEL.£E9.0) GU TD 33 
TKNC KK) we SMLE(L11) 
EXPM1 SE XPON(NMKM1{ 4,2) 
00 35 Qs1,NR80 
KTSKB4¢(Qet)aNRIFCH6 
KKSKT#NETTC 

TKN( KK) weTKNCKT) #E XPM 
DO 38 As1,NRBD 
KT3KB4(Qe1) eNRIFCOeNRAD 46 
KKK T#NETTC 

TKNCKK )BTKN(KT) wEXPMY 
GO TO 100 

Nisei 

N2@26 

DO 50 NN&@4,NE 
IF(NN,EQ,1)G0 TH 42 
N13] 

Nese 

NPMKSNPKONI 
INMKSBIABS(NPMK) 
RINMK8&1,0wINMK 

RF ABRINMK /NBP 

Nf& ABINMK /NBP 

DIF SABS(REA@1, ORNFA) 
IF COIF .G7,(,95))6U TN 50 
DO 44 Qai,NCOUOLS 
L4USLLARo (Qe )alae4u 
KKSKB+(()e] )RNRIFCHN2 
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44 


4§ 


46 


47 
50 


90 
100 


Liisl4¢7 

TKN(KK) @NARB(LO)ONI WCM1 eB(L11) 
L4aLSmaed 

KKSKBG4N2 

L1i1SL47 

IFCMCT,E0,0)G0 Th 45 

TRN( KK De eNARSMLE(LO)oNIMCMIWSMLB(L11) 
KKEKK#NRIFC 

IF CMFEA,EQ.0)GM TU 46 

TKN(KK )meNARSHLO(L4O)eNIMCMIMSMLC(L 11) 
KKEKK4#NRIFC 

IF(MFLAP,E9,0) GN TO 4? 

TKN(KK) SeNAWSMLO(L4) NI BCMIRSMLD(L11) 
KKSKK@NRIFC 

IF (MLEL,EQ,0)G0 TU SO 

TKN (KK) SeNARSMLE (LO) @NIRCMIERSMLE(L 11) 
CONTINUE 

GQ TO 100 

STOP 

RETURN 

END 


SuBRUUTINE POLAR 

CUMPLEX OXJ,QYI,QZJ,NTX,OTY,0TZ 
CUMMUN/SQVTEM/QOXJ oe OYJ, QZN 
DIXSIXTeCMPLX(0,0,01,0) 
DTYSOYJeCMPLX(0,0,01,0) 
DTZENZIeCMPLX(0,0,91,0) 
OXRBUXd 

DYRBOYJ 

NZRsOZJ 

DxleOTy 

OYISDTY 

DZI=0TZ 

IF (DXR,NE.0,0) GN TO 2 
TF (COXI,NE.0,0) GO TO 2 
DxAz0,0 

GU TO 3 

2 DXASATANZ(DXT,DOXR) 

3 IFCOYR,NE,0,0) GN TO 4 
IFCOYI,NE.0,0) GO TO 4 
DyYas0,0 
GU TO 5 

4 DYASATAN2(DYI,DYR) 

S IFCOZR,NE.0,.0) GO TO 6 
IF (OZI,NE,9,.0) GO TO 6 
D7A£0,0 
GO TO 7 

6 DZASATAN2(DZ1,0Z2R) 

7 CONTINUE 
DXRESORT(DXR#DXRODXT #DXIT) 
OYRESURT(CDYR®DYRODYI DVI) 
OZRASYRTCOZREDZROOZIaDZI) 
GXJeCMPLX(OXR, DXA) 
QYJSCMPLX(DYR,DYA) 
OZJBCMPLX(DZR,D0ZA) 

RETURN 
END 


*FORTRAN CALLABLE COMPLEX FUNCTION TN ORTAIN NAT PRONUCTS. 
*ARGUMENT LIST IS (NAB), WHERE N TS THE DIMENSTON OF THE VECTORS 
*A ANDO Be A IS PRESUMFD SPARSE FOR MAXIMUM PREGRAM SPFED. 
*INTE®MENTATE RFSULTS ARF CARRIFN IN NCUBLF PRECISION AND THE 
*FUNCTION MAY BF DECLAREN DOURLE PRECISION COMPLEX, IF OFSTRED. 

* 


SPACE 2 
#INCR EQU 0 
#COMPR EQU l 
#IANDEX EQU 2 
#N EQU 2 
#A EQU 3 
#R E9QU 4 
#MAXR FOU 4 
SPACE 
#REAL EQU 0 
#IMAG FOU 2 
#ZERQ EQU 4 
#TEMP EQU 6 
SPACE 2? 
A DSECT 
AREAL DS n 
AIMAG DS D 
B DSECT 
BREAL OS D 
BIMAG DS D 
EJECT 
cnor CSECT 


SAVE (2,#MAXR)D¢,* 
USING CNOT,15 
LM #N,#B,0(1) 
USING A,#A 
USING B,#B 
L #COMPR,O(AN) 
BCTR #COMPR,0O 
SLA #COMPR 94 
LA HINCR»16 
SR #INDEX »#INDEX 
SOR #REAL »#REAL 
SOR #IMAG,#IMAG 
SOR #ZERO,#ZERN 
SPACE 2 
LooP co #ZERO,A(#INDEX) 
BNE CONTINUE 
BXLE #INDEX,#INCR,LOOP 
EXIT RETURN (2 5#MAXR) 
SPACE 
CONTINUE LD #TEMP, AREAL ( #INDEX) 
*D #TEMP »-BREAL( #INDEX) 
ADR #REAL »#TEMP 


Ln #TFMP ,ALTMAG( 4TNDEX ) 
Mp) 4#TFMP,BIMAG(#INDEX) 
spe HREAL ,RTEMP 

LD #TEMP, AREAL ( #INNDEX ) 
Mn #TEMP ,BIMAG(#INNDEX) 
ANR #IMAG,4#TEMP 

Lo #TEMP,ATMAG(#INNFX) 
Mn #TEMP,BREAL (#INDEX) 


ADR 4#IMAG,#TEMP 

RXLF #INNDEX,#INCR Ll O0P 

8 EXIT 

SPACE 2 

LTORG 

ENN . ———-_ (Start of RCDOT) 
*FORTRAN CALLARLE COMPLEX FUNCTICN TO CRTAIN NOT PRONUCTS. 
*ARGUMENT LIST FS (NeAeB)», WHERF N TS THE DIMENSION OF THE VECTORS 
*A AND Be. A TS PRESUMED SPARSF FOR MAXIMUM PROGRAM SPEED. 
*A TS A REAL VECTOR, WHILE B IS COMPLEX. 
*INTERMENTATE RESULTS ARE CARRTFO IN NOQUBLF PRECISION AND THE 
*FUNCTION MAY PE DECLARED NOUBLF PRECISION COMPLEX, IF CFSIRED. 
*x 


SPACE 2 
#1 NCR EQU 0 
#COMPR EQU 1 
#INDEXA FOU 2 
#INNEXB EQU z) 
#N EQU 3 
# A FQU 4 
#8 EQU 5 
H#MAXR EQU 5 
SPACE 
#REAL EQU 0 
#IMAG EQU 2 
#ZERO FOU 4 
#TEMP EQU 6 
SPACE 2 
A NSECT 
B OSFCT 
BREAL OS D 
RIMAG oS D 
EJECT 


RCNNT CSECT 
SAVE (2,#MAXR) 9 e* 
USING RCNOT,15 
LM #N,#B,0(1) 
USING A,#A 
USING 8,48 
L #COMPR,OTHN) 
BCTR #COMPR,.O 
SLA #COMPR,3 


LA #INCR»8 

se #INTEXA,#INDEXA 
SOR HREAL,#REAL 

sor 4IMAG,HIMAG 


SOR #7ERN,#7ZFQN 
SPACF 2? 
LANe cod #ZEROA(#INDEXA) 


BNE CONTTNUF 
RYLE #INDEXA,#INCR,LENP 


FXIT PETURN (2,#MAXQ) 
SPACE 
CONTINUF LA HINDEXR,O(AINDFEXA,#INDEXA) | 
iD) H¥TEMDZALHTNDEXA) 
M1 #TEMP,BREAL (#INDEXR) 
ADR HREAL #TEMP 
LD #TEMP,ZA(HINODEXA) 
vy #TEMP »RIMAG(HINDEFYP ) 


ADR  #1MAG,#TEMP 

RXLE #INNEXA,#INCP,LONP 

R FXIT 

SPACE 2 

LTARG (Start of ROWSUM) 

END 2 7 

*FORTRAN CALLABLE SUBROUTINE TC PFREARM MATRIX @NW OPERATIONS. 
*ARGUMENT LIST IS (NeNOIM,A,R eX) 

*THE ROW NPFRATION A=A—X*R IS PFRFNORMFD, WHERE Ay By AND X ARE 
*NPUBLE PRECTSION COMPLEX, 

*NDIM IS THE COLUMA DIMENSTCN MF THE MATRICFS, AND N IS THE NUMBER OF 
*ELEMENTS TO RE OPFRATEN ON IN THF ROWS. THE INDEXING SCHEME IS 
*THEREFORE ACT) =ACTI-X*BU TD, T=1y LAEN-LI*NDIM,NOIM 

* 


SPACE 2 
#INCR FQU 0 
HCCMPR FOU 1 
#ITNDEX EQU 2 
4#N FOU 1 
ANN IM EQU 2 
#A EQU 3 
#9 FQU 4 
#X EQU 5 
AMAXR FOV 5 

SPACE 
#ATEMP EQU 0 
#RTEMP EQU 2 
#XRFEAL EFQU 4 
#X IMAG EQU 6 

SPACE 2 
A OSECT 
AREAL ns D 
AIMAG nS 0 


408 


a DSFCT 

BEFAL DS 0 

RIMAG MS nD 
FJECT 


ROWSUM = CSERT 
USING RAWSUM 415 
SAVE (14 #YAXRD Gg * 
LM 4 4#X,O(1) 
USING Ay#A 
USING B,#A 


L #NDIM,O(#NOTM) 
SLA 4#NNIM,4 
L #CAMPRAOCHN) 
RCOTR 4#COMPR,D 
MR #COMPR=-1,4NDIM 
LR #TNCR,#NDIM 
SR #INMEX,#TNNEX 
19 #XREAL,O( 4X) 
L9 #XIMAG,R( 4X) 
SPACE 2 
Lare ia) HATEMD,ARIMAG(HINNEX) 
MOR H#ATEMP HX IMAG 
AND #ATEMP,ARFAL (#TNDEX) 
19 H#ATEVP,RREALTH#INDFEY) 


MNR #3TEMP,HKPFAL 
SOR #ATEMP,H#BRTEMP 
STO HATEMP AREAL (#INNFX ) 
TD) #ATEMP ,ATMAGIHINDFX) 
ap) #RTEMP,BREAL (HINNFY) 
MIR H8TEMP,#XTMAG 
SNR #ATEMP,HRTEMP 
aD) H#ATEMP,PTMAGIHINNEYX 
MPR HATEMP ,#XRFAL 
SNR HATEMD,#R TEMP 
stn HATEMP,AIMAG(H#INDFEX) 
BXLE #INDEX,#INCR,LOANP 
SPACF 
RETURN (2 9#MAXR) oT 
LTORG 
FENN 
OVERLAY ALPHA (Start of Overlay Structure) 
INSERT ART 
INSERT ARR 
INSERT COEFFS 
INSERT BTS 
INSFRT SS] 
INSEPT TKNI 
OVERLAY BETA 
INSERT BAFRO 
INSERT FAERO 
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INSERT AFRO 
INSERT ACIEFF 
INSERT SETUP 
INSERT SFCPAQ 
INSFPT 1 AAOIN 
INSERT TAALU 
MOVERIAY RETA 
INSFAT FLAST 
INSEQT RIGIN 
INSFOT STIFF 
INSFOT BEND 
INSFRT MLeC2 
INSERT MLCC2 
INSF®8T CHOT 
INSEOT 2COAT 
AVERLAY GAMMA 
INSERT TARRAY 
INSERT RASS 
INSERT BEARN 
NOVFRLAY GAMMA 
JINSFRT SAQDAY 
INSE2T EMASS 
INSFET FUARN 
AVERLAY PETA 
INSFOT FROSOEN 
INSE&T SUBMA 
INSFPET SURMR 
INSERT SUBME 
INSERT SUBYG 
INSFRT 7TEGT 
INSF2T FXPON 
INSERT TENS 
INSFRT SWwA 
INSFAT Swe 
INSERT 7IN 
INSFRT X¥NEO 
INSFEPT ULN 
INSFPT FYCHT 
OVERLAY ALPHA 
INSFOT SALVE 
INSERT OCMAT 
TNST R&T SWAPS 
INSE2T ROWSUM 
INSFaT FRRSET 
ENTOY MAIN 


aie ee ee 
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